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ABSTRACT 

Accurate  modelling  of  gun  interior  ballistics  promotes  more  efficient  gun  and 
propelling  charge  design.  In  order  to  simulate  interior  ballistic  flowfields, 
such  models  require  a  description  of  the  thermodynamic  behaviour  of  the 
propellant  gas.  The  Noble-Abel  equation  provides  a  simple  and  reasonably 
accurate  equation  of  state  for  propellant  gases  at  the  high  densities  and  tem¬ 
peratures  experienced  in  guns.  Most  computational  fluid  dynamics-based  bal¬ 
listics  models,  however,  require  additional  thermodynamic  functions  which 
must  be  derived  from  the  equation  of  state.  This  note  presents  the  derivation 
of  such  thermodynamic  functions  for  Noble-Abel  gases.  Although  the  deriva¬ 
tions  are  geared  toward  the  functional  requirements  of  the  commercial  Fluent 
code,  the  results  are  equally  applicable  to  all  computational  fluid  dynamics 
solvers.  Also  presented  is  a  brief  numerical  example  for  a  typical  propellant, 
highlighting  the  different  thermodynamics  of  the  Noble-Abel  and  ideal  gas 
equations. 
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The  Noble-Abel  Equation  of  State:  Thermodynamic  Derivations 

for  Ballistics  Modelling 


EXECUTIVE  SUMMARY 

Accurate  modelling  of  gun  interior  ballistics  promotes  more  efficient  gun  and  pro¬ 
pelling  charge  design.  In  order  to  simulate  interior  ballistic  flowfields,  such  models  re¬ 
quire  a  description  of  the  thermodynamic  behaviour  of  the  propellant  gas.  The  Noble- 
Abel  equation  provides  a  simple  and  reasonably  accurate  equation-of-state  for  propellant 
gases  at  the  high  densities  and  temperatures  experienced  in  guns. 

Most  computational  fluid  dynamics-based  ballistics  models,  however,  require  addi¬ 
tional  thermodynamic  functions  which  must  be  derived  from  the  equation-of-state.  This 
note  presents  the  derivation  of  a  range  of  thermodynamic  functions  for  Noble-Abel  gases. 
They  include: 

•  Entropy, 

•  Speed  of  sound, 

•  The  functional  form  of  the  specific  heats, 

•  The  relationship  between  the  specific  heats  and  gas  constant, 

•  The  isentropic  process,  and 

•  Various  partial  derivatives  of  density  and  enthalpy. 

Although  the  derivations  are  geared  toward  the  functional  requirements  of  the  commer¬ 
cial  Fluent  code,  the  results  are  equally  applicable  to  all  computational  fluid  dynamics 
solvers. 

Also  presented  is  a  brief  numerical  example  for  a  high-energy  tank  gun  propellant 
(JA2),  demonstrating  the  difference  between  the  Noble-Abel  and  ideal  gas  equations-of- 
state  and  two  of  the  derived  functions. 
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Nomenclature 

Speed  of  Sound  [m/s] 

Co-volume  [m3/kg] 

Specific  Heat  at  Constant  Pressure  [J/(kg  K)] 

Specific  Heat  at  Constant  Volume  [J/(kg  K)] 
Intensive  Total  Energy  [J/kg] 

Internal  Energy  [J/kg] 

Vector  of  Fluxes 
Enthalpy  [J/kg] 

Normal  Unit  Vector 
Pressure  [Pa] 

Specific  Gas  Constant  [J/(kg  K)] 

Entropy  [J/(kg  K)] 

Surface  Area  [m2] 

Time  [s] 

Temperature  [K] 

Velocity  Vector  [m/s] 

Vector  of  Conserved  Variables 
Specific  Volume  [m3/kg] 

Volume  [m3] 

Coefficient  for  Intermolecular  Attraction  [m5/(kg  s2)] 
Ratio  of  Specific  Heats 
Density  [kg/m3] 
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1  Introduction 


The  capability  to  accurately  model  gun  interior  ballistics  (IB)  promotes  faster  and 
cheaper  gun  and  propelling  charge  design.  Likewise,  modelling  can  be  used  to  predict 
and  assess  a  gun's  performance  characteristics  before  it  is  built,  tested,  purchased  or  mod¬ 
ified.  Although  such  a  predictive  utility  can  be  partially  achieved  through  targeted  ex¬ 
periments,  models  provide  a  range  of  additional  benefits.  The  extreme  interior  ballistic 
environment  makes  experimental  instrumentation  and  measurement  difficult,  whereas  it 
is  normally  possible  to  determine  all  modelled  physical  quantities  throughout  a  simula¬ 
tion  domain.  The  ability  to  add  or  remove  different  physical  phenomena  at  will,  allows 
models  to  be  used  to  identify  the  relative  action  of  the  various  ballistic  processes.  Mod¬ 
elling  also  allows  automated  optimization.  For  example,  charge  weight  can  be  minimized 
by  optimizing  propellant  grain  geometry  subject  to  muzzle  velocity  and  maximum  pres¬ 
sure  constraints. 


One  of  the  simplest  classes  of  IB  models  are  lumped  parameter  models,  such  as  IB- 
HVG2  [1],  where  the  dynamic  firing  process  is  represented  by  mean  (lumped)  state  vari¬ 
ables.  A  particular  pressure  gradient  between  breech  and  projectile  base  is  assumed 
and  applied  to  calculate  projectile  base  pressure  and  its  resulting  acceleration.  Two- 
phase,  multi-dimensional  computational  fluid  dynamics  flow  (CFD)  solvers,  such  as  the 
NGEN  [bj]  code,  provide  higher-fidelity  IB  simulations  but  are  computationally  expen¬ 
sive.  Physical  phenomena  such  as  inter-phase  drag,  axial  and  radial  flame  spreading 
within  the  propellant  bed,  boundary  layer  formation,  and  pressure  waves  can  be  mod¬ 
elled  using  this  class  of  tools.  Of  intermediate  complexity  are  codes  such  as  XKTC  [j3[|, 
which  combines  a  one-dimensional  flow  solver  with  a  lumped-parameter  model. 

Irrespective  of  which  approach  is  used,  all  models  require  an  accurate  description  of 
the  thermodynamics  of  the  propellant  gas.  Consider  the  Euler  equations,  which  may  be 
thought  of  as  a  simplified  version  of  the  Reynolds-averaged  Navier-Stokes  equations  that 
are  used  by  the  CFD  solvers  to  describe  the  ballistic  flow.  Assuming  no  source  terms,  for 
a  control  volume  V  the  Euler  equations  can  be  written  in  integral  form  as 


i,Jviv+ 


=  0, 


(1) 


where  U  is  a  vector  representing  conserved  flow  quantities  at  points  within  the  control 
volume,  F  is  a  vector  of  fluxes  across  the  surface  of  the  control  volume,  the  scalars  V  and 
S  refer  to  the  volume  and  surface  of  the  control  volume  respectively,  and  t  is  the  time 
variable.  For  a  single-species,  single-phase  flow,  the  vectors  U  and  F  may  be  expressed  as 


p 

P( u  •  n) 

u  = 

pu 

,  and  F  = 

pu(u  •  n)  +  Pn 

IPE] 

pE( u  •  n)  +  P(u  •  n) 

Here,  u  is  the  fluid  velocity  vector  and  n  is  a  unit  normal  to  the  control  volume  surface.  The 
primitive  variables  p,  P,  and  E  represent  density,  absolute  static  pressure,  and  intensive 
total  energy  respectively.  Total  energy  can  be  expressed  in  terms  of  internal  energy  and 
kinetic  energy  as 

E  =  e+\  |u|2,  (3) 
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and  internal  energy  can  be  calculated  by  integrating  the  constant-volume  specific  heat 
with  respect  to  temperature: 


(4) 


Equations  [g-  a\  contain  six  variables  but  provide  only  five  relations.  To  solve  the  system 
an  additional  equation— relating  the  state  variables  of  the  fluid— must  be  provided. 


For  gases  at  moderate  to  low  density,  the  ideal  gas  equation  of  state 


P  =  pRT, 


(5) 


where  R  is  the  specific  gas  constant,  can  be  used  to  close  the  system.  The  ideal  gas  equa¬ 
tion  of  state  is  accurate  providing  that  (1)  the  average  intermolecular  spacing  or  mean 
free  path  of  the  gas  is  very  large  compared  to  the  size  of  the  gas  molecules,  and  (2)  inter¬ 
molecular  forces  are  weak.  For  many  practical  flow  regimes  both  of  these  assumptions 
are  true.  However,  the  high  gas  densities  occurring  in  a  gun  chamber  during  combustion 
render  the  ideal  gas  equation  of  state  inaccurate.  Solid  loading  densities  for  gun  pro¬ 
pelling  charges  are  typically  of  the  order  500-1000  kg/m3.  Peak  gas  densities  produced 
during  the  ballistic  cycle  may  also  approach  this  magnitude  since,  in  most  guns,  the  bulk 
of  propellant  is  burnt  before  the  projectile  has  moved  very  far  downbore. 

While  still  an  approximation,  the  van  der  Waals  equation  of  state  provides  an  im¬ 
provement  in  accuracy  for  high  density  gases.  It  is  of  the  form  [I4J] 

(P  +  u/v2){v-b)  =  RT,  (6) 


where  the  gas  specific  volume  v  =  1/p.  The  co- volume  b  compensates  for  the  finite  vol¬ 
ume  occupied  by  the  gas  molecules,  while  the  term  oc/v 2  accounts  for  intermolecular  at¬ 
traction  forces.  Note  that  in  the  case  a  =  b  =  0,  Equation  6  reverts  to  the  ideal  gas  equation 
of  state.  In  ballistics  applications,  the  high  propellant  gas  temperature  means  that  inter¬ 
molecular  attraction  energy  is  small  in  comparison  to  molecular  kinetic  energy  II4J].  Thus 
the  attraction  term  can  be  removed  without  significant  loss  of  accuracy,  resulting  in  the 
so-called  Noble-Abel  equation  of  state: 

P(v  -b)  =  RT.  (7) 


Equation  ^]can  be  used  to  describe  the  propellant  gas  with  sufficient  accuracy  for  both 
lumped  parameter  and  CFD  IB  models.  For  lumped  parameter  modelling,  nothing  fur¬ 
ther  is  required.  Most  CFD  models,  however,  require  a  range  of  additional  thermody¬ 
namic  functions  corresponding  to  this  equation  of  state.  It  is  necessary  to  calculate  sound 
speed,  for  example,  in  order  to  evaluate  the  Courant-Friedrich-Fewy  criterion  and  main¬ 
tain  numerical  stability;  the  expression  for  speed  of  sound  in  an  ideal  gas  is  different 
to  that  for  Noble-Abel  gases.  Similarly,  partial  derivatives  of  state  variables  may  be  re¬ 
quired  quantities  for  certain  implicit  CFD  schemes.  The  commercial  Fluent  CFD  code  I5I], 
for  example,  requires  eight  thermodynamic  functions  that  must  be  derived  from  (or  are 
related  to)  the  equation  of  state.  The  aim  of  this  note  is  to  present  the  derivation  of  such 
thermodynamic  functions  for  a  Noble-Abel  gas.  Although  this  note  closely  follows  the 
requirements  of  the  Fluent  code,  the  results  are  equally  applicable  to  all  commercial  and 
research  CFD  solvers.  Also  presented  is  a  numerical  example  for  a  typical  propellant, 
comparing  thermodynamic  quantities  calculated  under  the  Noble-Abel  equation  of  state 
with  those  of  an  ideal  gas. 
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2  Thermodynamic  Functions 


The  Fluent  code  requires  definitions  for  the  thermodynamic  functions  listed  below  Jig]. 
These  functions,  and  some  additional  relations,  are  derived  in  Sections  2. 1-2.7. 

p(T,P)  Cp(T,p)  h(T,p )  s(T,p)  a(T,p ) 


dp 


dp 


dP 


i f)(T'p )  (^)om  (sfJM 


dh 


dT 


2.1  Equation  of  State 

The  Noble-Abel  equation  of  state  has  already  been  defined  in  the  introduction  as 

P(v  -b)=  RT, 

however  Fluent  requires  an  equation  of  state  in  the  form  p(T,P).  It  is  easy  to  recast  the 
equation  as 

p  =  P/{RT  +  bP )  (8) 

to  satisfy  this  requirement. 


2.2  Specific  Heats 


Fluent  also  requires  the  specific  heat  at  constant  pressure  as  a  function  of  density  and 
temperature.  For  completeness  we  will  consider  both  specific  heats  in  this  subsection. 
They  are  defined  as 


(  dh\ 


fde\ 


Cv  =  ^  and  cz,  =  — 


\dT  J 


\dT  J , 


(9) 


From  Reference  [g],  the  gradient  of  Cp  with  respect  to  pressure  at  constant  temperature 
can  be  related  to  the  equation  of  state  via 


d  cp 
~dP 


=  —T  I 


(  d2V  \ 


(10) 


Evaluating  the  right  hand  side  for  the  Noble-Abel  equation  of  state  yields 

dcv\  m  d2  (RT  ,\ 

ar)T  =  -T3pU  +  i,)P  =  a  (II) 

By  Equation  [8]  this  also  proves  that  (dCp/dp)T  =  0,  and  hence 

cp{P'T)  =  cp(T)  (12) 

for  a  Noble-Abel  gas.  A  similar  proof  is  possible  for  specific  heat  at  constant  volume: 

RT 


dc,\  = 

dv)T  \d  T2)r,  d  T2\v-b 


=  0, 


(13) 
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and  thus 

cv(p,T)  =  cv(T).  (14) 

Results  12  and  14  show  that,  for  a  Noble-Abel  gas,  the  specific  heat  functions  can 
be  conveniently  implemented  in  computer  code  as  a  curve  fit  or  look-up  table  in  terms 
of  temperature  only.  As  an  alternative  to  providing  the  specific  heats  as  a  function  of 
temperature,  the  calorifically  perfect  approximation  could  also  be  made: 

Cp  =  constant  and  cv  =  constant.  (15) 

2.3  The  Specific  Gas  Constant 

The  previous  section  showed  that  both  c;,  and  cv  could  be  expressed  as  functions  of 
temperature  only.  We  now  proceed  to  show  that  these  quantities  can  be  related,  such  that 
if  one  specific  heat  is  known  then  the  other  may  be  calculated  easily.  This  obviates  the 
need  to  provide  two  separate  tables  (or  curve  fits).  From  Reference  [J6J 


Evaluating  the  partial  derivatives  from  the  Noble-Abel  equation  of  state  yields 

_t{R\2  RT 

c”  C'-T(pJ  (v-br 

Upon  simplifying,  we  get  the  convenient  relationship 

Cp  Cd  —  R, 

which  is  the  same  as  that  for  ideal  gases. 


(17) 


(18) 


2.4  Entropy 

Entropy  for  a  Noble-Abel  gas  may  be  derived  starting  from  the  thermodynamic  rela¬ 
tion 

Tds  =  de  +  PdV.  (19) 

Substituting  with  Equations  j7]  and  [9}  we  have 

dTR 

ds  =  cv—  - 1 - -dv.  (20) 

T  v  -  b 

Equation [20] may  be  integrated  to  calculate  the  entropy  difference  between  two  states: 

s2  -  si  =  r  c^dT  +  r  — ~dv .  (21) 

Jt !  T  JVl  v-b 
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The  second  term  is  easy  to  evaluate,  however  the  first  term  can  only  be  analytically  inte¬ 
grated  if  the  form  of  cv  is  known.  For  the  case  of  a  calorifically  perfect  gas  (Equation  15), 
integration  yields 

or  equivalently,  in  the  form  required  by  Fluent, 

S2-S1  =  c„ln  (5)  +Rln  •  te) 


2.5  The  Isentropic  Process 


Although  not  directly  important  to  ballistics  modelling,  the  consideration  of  state 
changes  occurring  at  constant  entropy  yields  a  result  that  will  be  of  use  in  the  following 
section.  For  an  isentropic  process  Si  =  S2,  and  for  a  calorifically  perfect  gas  Equation  22 
gives 


and  thus 


Ti 

Ti 


--In 

f  v2  -  b\ 

(24) 

V^i  -  b) 

Cv 

Vi  -  b) 

R/cv 

(25) 

v2  -  b  ) 

1 

Replacing  the  temperature  terms  by  pressure  and  specific  volume  (Equation  [7J  yields 


Pi(v  2-b)  =  Ah-^\R/Ci 
Pi(vi-b)  \v2-bj 


(26) 


and 

Pi {vi  -b) ( 1+R/C =  P2(v2-b) (1+R/c") .  (27) 

The  ratio  of  specific  heats  is  defined  as  7  =  cp/cv,  which  together  with  Equation  18  results 
in  the  expression 

P(v  —  b)7  =  constant  (28) 

for  an  isentropic  process  in  a  calorifically  perfect,  Noble-Abel  gas.  This  is  analagous  to 
the  familiar  ideal-gas  expression  Pv 7  =  constant. 


2.6  Sound  Speed 


The  speed  of  sound  in  a  gas,  a,  is  defined  as  the  speed  of  reversible  pressure  wave 
propagation.  From  Reference  p],  this  is 


a  = 


dP\ 


1/2 


We  may  expand  the  partial  derivative  using  the  chain  rule 


dp\  _  /dPdv\ 


(29) 


(30) 
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and  since  specific  volume  is  just  the  reciprocal  of  density 


(31) 


We  now  take  the  partial  derivative  of  both  sides  of  Equation^  with  respect  to  specific 
volume 

^[p(v-by]  =  0,  (32) 

and  by  the  product  rule 


—  (v  —  b)7  +  Py(v  —  b)^7  =  0. 

av 


(33) 


Remembering  that  Equation  [28] represented  an  isentropic  process  in  a  calorifically  perfect 
gas,  we  therefore  have 


dP\  _  P  7 


dv  J , 


v  —  b 


(34) 


Finally,  substitution  of  Equations  [34] and  [7] into  Equation  [31] yields  the  speed  of  sound  in 
a  calorifically  perfect  Noble-Abel  gas: 


a  =  — V—r^yRT. 
v  -  b  v 


(33) 


This  expression  may  be  thought  of  as  the  speed  of  sound  in  an  ideal  gas,  y ^yRT,  with  a 
premultiplying  coefficient  to  correct  for  the  effect  of  co-volume.  It  is  easily  rewritten  in 
the  form  required  by  Fluent, 


(36) 


It  can  be  shown  that  Equation  ^36]  also  holds  for  the  more  general  case  of  Noble-Abel 
gases  which  are  not  calorifically  perfect.  In  Reference  Anderson  derives  the  general 
expression 


a 


2 


7 P  ri  +  (l /P)(de/dv)T 
~p  [  l-p(dh/dP)T 


(37) 


for  the  sound  speed  of  non-ideal  gases,  reacting  gases,  and  thermally  and  calorifically 
imperfect  gases.  The  partial  derivatives  of  Equation  [37]  can  be  evaluated  using  two  ther¬ 
modynamic  relations  given  in  []6[|.  The  partial  derivative  of  internal  energy  with  respect 
to  specific  volume  at  constant  temperature  can  be  related  to  the  equation  of  state  via 


^  =T(d-l)  -P. 
dv)T  \dT  J  v 


(38) 


Upon  substition  of  Equation  [7}  we  have 


/de' 

[dVy  T 


RT 
v  —  b 


-P  =  0. 


(39) 
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Likewise,  the  partial  derivative  of  enthalpy  with  respect  to  pressure  at  constant  tempera¬ 
ture  can  be  evaluated  by  the  relation 


\dP  J  T 

For  the  Noble-Abel  equation  of  state,  this  yields 

( dh  \  RT  , 

(arjT  =  !’- T  =  ‘- 

After  substitution  of  Equations  gi]  and  39  into  Equation  [37]  we  have 

1 


„2=7P 


which  can  be  simplified  to  give 


a  = 


p  [l  — bp 


—P^RT- 


(40) 


(41) 


(42) 


(43) 


This  expression  is  identical  to  the  sound  speed  for  calorifically  perfect  Noble-Abel  gases 
given  by  Equation  36 


2.7  Partial  Derivatives  of  State  Variables 


The  Fluent  code  additionally  requires  that  various  partial  derivatives  of  state  variables 
be  supplied.  Derivations  of  the  non-trivial  partial  derivatives  are  now  provided. 

First  we  the  derive  the  partial  derivative  of  density  with  respect  to  pressure  at  constant 
temperature,  starting  with  the  chain  rule 


fdp)  =  (dpdv\ 

\dPjj  \dvdP  J  T' 

and  thus 

(dp\  =_L(dl\ 
\dP  J  T  v2\dPjT 


From  the  Noble-Abel  equation  of  state 


so 


-  b)T 


fdv\  _  RT 
{ dp)T~~P I' 

Substituting  this  into  Equation  45  and  simplifying  yields  the  result 


(44) 


(45) 


(46) 

(47) 


(48) 
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which  can  be  rewritten  in  the  required  form  as 

^  (1  -pb)2 
T 


dP 


RT 


(49) 


Fluent  also  needs  the  partial  derivative  of  density  with  respect  to  temperature  at  con 
stant  pressure.  A  derivation  similar  to  the  above  gives 

dp  \  _  bp2  —  p 

dT  L  =  T  ' 


(50) 


Finally,  two  partial  derivatives  of  enthalpy  are  required  by  Fluent.  The  first,  ( dh/dT)p 
is  by  definition  equal  to  cp  (Equation  9).  The  second,  ( dh/dP)j  is  equal  to  the  co-volume 
b  (Equation  fci). 


3  Example 

A  brief  example  is  now  presented,  to  demonstrate  the  difference  between  the  Noble- 
Abel  and  ideal  gas  equations  of  state  and  two  of  the  derived  functions.  A  high  perfor¬ 
mance  tank  gun  propellant,  JA2,  will  be  used  as  the  subject.  The  approximate  thermody¬ 
namic  properties  of  JA2  propellant  gas,  including  its  specific  heats  at  the  propellant  flame 
temperature,  are  shown  in  Table 

Figure  2]  shows  pressure  as  a  function  of  density  and  temperature  for  JA2  propel¬ 
lant  gas,  according  to  the  Noble-Abel  equation  of  state.  The  variation  between  pressure 
calculated  using  the  Noble-Abel  and  ideal  gas  equations  is  represented  by  the  coloured 
contours.  At  a  typical  IB  gas  density  of  300  kg/m3,  for  example,  the  pressure  difference 
is  30%.  The  differences  in  pressure  are  observed  to  be  independent  of  temperature.  This 
is  to  be  expected,  since  the  ideal  and  Noble-Abel  equations  of  state  can  be  related  by  a 
substitution  in  terms  of  specific  volume  (or  density)  only:  v  <->  (v  —  b). 

Figure  [2]  shows  sound  speed  as  a  function  of  density  and  temperature  for  JA2  pro¬ 
pellant  gas,  assuming  that  it  is  calorifically  perfect  and  obeys  the  Noble-Abel  equation  of 
state  (Equation  35).  Again,  the  difference  between  the  Noble-Abel  and  ideal-gas  values 
are  represented  by  coloured  contours  and  depend  on  density  only.  The  differences  are  of 
similar  magnitude  to  those  already  observed  for  pressure. 

Entropy,  as  a  function  of  density  and  temperature,  is  shown  in  Figure  g]for  the  Noble- 
Abel  equation  of  state.  The  gas  is  assumed  calorifically  perfect.  A  reference  entropy  (si  in 
Equation [22)  corresponding  to  a  state  of  300  K  and  1  kg/m3  has  been  used.  The  difference 
between  the  Noble-Abel  and  ideal-gas  entropy  is  shown  in  the  figure  in  terms  of  absolute 
entropy  values. 

The  effect  of  the  equation  of  state  on  predicted  ballistic  performance  is  dependent 
on  the  particular  gun  and  propellant  properties,  and  the  operating  regime.  In  the  case  of 
JA2  propellant,  muzzle  velocity  predictions  for  typical  high  performance  medium-calibre 
guns  can  be  expected  to  differ  by  around  10-15%.  Likewise,  peak  gun  chamber  pressure 
predictions  will  be  significantly  higher  (around  30-35%)  when  the  Noble-Abel  equation 
is  used  instead  of  the  ideal  gas  equation.  This  is  consistent  with  the  results  of  Figure  ,T 
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Table  1:  Properties  ofJA2  Propellant  Gas 


Gas  Constant,  R 

334  J/(kg  K) 

Specific  Heat  Ratio,  7 

1.225 

Specific  Heat  at  Constant  Volume,  cv{Tf) 

1484  J/ (kg  K) 

Specific  Heat  at  Constant  Pressure,  cp(Tf) 

1818  J/(kg  K) 

Co-volume,  b 

0.001  m3/kg 

Flame  Temperature  of  Solid  Propellant,  Tf 

3410  K 

Difference  [%] 

I  46 
42 
38 
34 
30 
26 
22 
18 
14 
10 
6 


Figure  i:  Pressure  as  a  function  of  density  and  temperature  for  JAi  propellant  gas,  according  to 
the  Noble-Abel  equation  of  state.  Coloured  contours  indicate  the  difference  between  ideal  gas  and 
Noble- Abel  gas  residts. 
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Figure  2:  Sound  speed  as  a  function  of  density  and  temperature  for  JA2  propellant  gas,  according 
to  the  Noble- Abel  equation  of  state.  Coloured  contours  indicate  the  difference  between  ideal  gas 
and  Noble-Abel  gas  results. 


Figure  3:  Entropy  as  a  function  of  density  and  temperature  for  JA2  propellant  gas,  according  to 
the  Noble-Abel  equation  of  state.  Coloured  contours  indicate  the  difference  between  ideal  gas  and 
Noble-Abel  gas  calculations. 
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